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Abstract 

Hadron decay widths are shown to increase in strong magnetic 
fields as F(eB) ~ ^T(O). The same mechanism is shown to be present 
in the production of the sea quark pair inside the confining string, 
which decreases the string tension with the growing eB parallel to the 
string . On the other hand, the average energy of the qq holes in the 
string world sheet increases, when the direction of B is perpendicular 
to the sheet. These two effects stipulate the spectacular picture of the 
B dependent confinement and a Sl discovered on the lattice. 
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The QCD confinement (as well as perturbative gluon exchange) was shown to 
be created by the nonperturbative (np) color-electric field correlators [1, 2] [3] 
which are not affected by magnetic field (m.f.) in the lowest order in a s . 
However, in the next order in a s (or in the 1/N C expansion) both confinement 
and gluon exchange (GE) interaction contain quark loops, which interact with 
the m.f. and can influence the resulting potentials. 

For the GE part it was found in [J], that the energy growth of the virtual 
qq in m.f. prevents the original QQ system from the collapse, keeping the 
GE interaction (VgeM)) finite at all eB. 

An interesting picture has emerged from the recent lattice studies in [5], 
where it was shown, that confinement interaction V con {(R) decreases for B 
parallel to R and increases for the perpendicular orientation, while \Vqe\ 
behaves in the opposite way. In the present paper we suggest an explana¬ 
tion of the behavior V con f(R ) and a s , and simultaneously we point out the 
stimulating role of m.f. in the strong hadron decay process. 
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The paper is organized as follows. In the next section we display the 
path integral Hamiltonian in m.f., the resulting wave functions, and some 
properties of the spectrum for the opposite charge qq systems. In section 3 
we derive shortly the magnetic focusing effect in the creation of the qq pair. 
In section 4 we describe the appearance of sea quark holes in the confining 
film and the m.f. dependence of the resulting effective string tension. We 
also discuss the dependence on the relative direction of m.f. and make a 
comparison with lattice data. In section 5 the a s dependence on m.f. is 
derived and compared to lattice data. In section 6 we compare our results 
with the effective action expansion and lattice data on average field strength 
squared in m.f. In section 7 a quantitative comparison of our results with 
the lattice data is presented. Section 8 is devoted to the summary of results 
and possible developments of the effects presented in the paper. 


2 Hamiltonian technic for hadrons in mag¬ 
netic field 


In this section we exploit the path-integral Hamiltonian approach for the 
qq systems (mesons) in m.f., which for the neutral case is embodied in the 
Hamiltonian O El IS] 


H 


p2 

2(u>i + U 2 ) 



+ U(v) + E 

i=l,2 


+ Uj 
2uji 


( 1 ) 


where P is the total momentum, n — and 77 is the relative qq distance, 
while mi — ^2 = m is the current quark mass. Here cUj is the (virtual) 
energy of the quark i which should be found from the minimum of the total 
energy eigenvalue, H\ J/ n = E n 'V n and 


dE n {u> u u> 2 ) 


duji 


,( 0 ) 


0 . 


( 2 ) 


The resulting stationary value E n (o;j; 0 \ ) is the actual energy of the qq 

system. 

I 11 the course of the decay the qq pair appears nearby the string connecting 
the original quarks Q and Q, which are assumed to be heavy for simplicity. 
We shall show in this section, that the magnetic focusing effect [2j is acting 
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since both q and q are charged. In this case (ignoring the c.m. motion, 
P = 0) one should write the total Hamiltonian as 


rr(<? 9 ) 

n QQ 


Pi + m{ 
2cai 


pi + ™% 

2u 2 


H---h Uiji — Rq) + U( r 2 — Rq). 


(3) 


The solution is readily obtains as a sum of two heavy-light mesons, centered 
at Rq and Rq. However, one should impose the condition of the relative 
state quantum number for qq which can be created by the nonperturbative 
(n.p.) or perturbative mechanism, yielding J PC = 0 ++ ( 3 Po mechanism) or 
1 ( 3 S'i mechanism) qq states respectively. 

Now let us switch on the m.f. The Hamiltonian ([3]) transforms as follows 


h qq(B) = E 


IP; 


e*A) 2 + oj 2 + m 2 


e t c t 


(i) B +(pl 


i=q,q 


2'jjj 


-U(t 1 -R q )+U(t 2 -Rq). 


(4) 

Here A(r) = |(B x r), and it is convenient to choose the origin r = 0 just in 
the middle of the distance (Rq — Rq). Now one can separate out the center 
of mass motion using the coordinates 


P 


coqri + uj 2 r 2 

LUl + u 2 


rj = r 1 - r 2 , 


ld_ 

i dr] ’ 


P 


ld_ 

i dp ’ 


(5) 


and one has 

H { q %\b) = H(P) + H w + U (6) 

where H( P) can be eliminated using the pseudomomentum procedure as in 
[8], U stands for the last two terms in (j3J), and 


= 


7r 

2 u 


E 

i= 1,2 


m, 2 + m 2 


e,cr,B 


2 a.’, 


+ 


CjViir] 2 


x 




(7) 


where fTO], (ei — e — — e 2 ) and subscripts _L and || refer to the direction of 


B, 


X v = 


e(uj 2 -Wi) 


n, = ^. 

* 2u 


2(jJ\uj 2 

We take into account, that expanding U in powers of one has U = <jR 
arjn + 0{rf) where R = Rq — Rq, n = and therefore disregarding U 
the first approximation, one has a solution for H n (uji = u> 2 = u>) 


( 8 ) 


M(u) = 


m 2 + uj 2 

UJ 


eB 

H-( 2n± + 1), 

UJ 


(9) 
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UJ=UJQ 


yielding at the stationary point ^ 


= 0, 


M 0 = M(ujq ) = 2^/m 2 + e q B. (10) 

Note, that in both cases 3 Po(S = 1 ,L = 1) and 3 S'i the spin and orbital 
projections cancel. Hence the qq pair acquires the effective mass (ITU , which 
grows with eB, when the qq loop stays in the confining him, when B is 
perpendicular to R. 

One can easily see in (17)) . (18)) that the situation is different in the case 
when B is parallel to the qq loop trajectory since in this case = 0 in ()8l) 
and the resulting M 0 = 2m in (fTOl) . 

However for the transverse m.f. the string acquires additional energy M 0 , 
Eq. dlOD, and the total energy of the QQ string with the qq hole can be 
estimated as E(R) = Vqq(R) + M 0 = oR + M 0 , and the resulting ratio of 
the energy increase per one hole is 


A E(B) _ 2^m 2 q + e q B 
E(R) = ~aR 


( 11 ) 


3 Magnetic focusing in the qq pair creation 


Magnetic focusing was treated in |9] in the case of two elementary objects; 
we now take the case of hadron constituents in ()4j) . Consider the expansion 
of U in the powers at the ratios i = 1, 2. Taking p = ri ^ r2 , one has 

U = <tR- <tv\\ + ^( 2 ^i + ^ 2 ), ( 12 ) 

where the subscripts (||) and (_L) stand for parallel and perpendicular with 
respect to R. Taking into account ()TJ) , and solving Htp = (H w + U)tp, one 
obtains the R-dependent wave function (for n± = 0) 


Mv±) 


exp(—r]'j_/2rY) 1 

\Eir ± ’ r]_ 



+ ( e q B) 2 , 


(13) 


where Uq is to be found from the stationary point of the ^-dependent energy, 
as in ()2I). For the lowest energy state one has from (1131) and (J7)). 


2cr 1 2au> 


E(oS) =ou + -\ — H— 


2 V coR 


00 


(e q B) 2 . 
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R 


(14) 














Taking the minimum of (ITTjh one finds r\ and hence </> 0 (0). Now the qq 
pair creation is described by the qq Green’s function [6] (in the background of 
the original QQ Wilson loop), G qq (x,y) ~ (x, x|e^^^ 4_y4 ^|yy) ~ \cp(0)\ 2 e~ ET . 
Hence the change in the wave function due to B can be characterized by the 
magnetic focusing factor 


, •pl(m = 0;eB) + ( e , B )2 

s . .2/ 


Ao(v± = 0; 0) 
One can see in (IT5l) two limiting cases 

a) e q B <C k 2 , £ = 1 + 


2cra;o(0) 

R 


(e,Bf 


R 2 = 


'\/3< 


a 


2/3 


R 


R 


(15) 


(16) 


b ) e q B > r 2 , £ ~ 6qB 


R 


2 • 


(17) 


For R « lfm one has r 2 = 0.14 GeV 2 (r 2 = 0.22 Gev 2 for R = 0.5 fm), 
and one obtains a strong amplifying factor for B||R and e£> ~ 1 GeV 2 . 


4 Sea quark effects in the confinement regime 


It is clear, that m.f. acts on the fixed boundary Wilson loop Wqq(A) through 
the creation of sea quark loops, which effectively create the holes in the him, 
covering the original Wilson loop. 

Following PHH2! one can write the partition function with the account 
of sea quark loops as 


Z = y DAexjp CAWQQ(A)det(m q + D(A). 


(18) 


where det{m q + D(A) can be written in the path integral form 

1 f°° ds 


det(m q + D(A )) = exp 


tr 


2 Jo s 


(D 4 z) e~ K ^W qq (A) 


Here W qq is the closed loop of the sea quark and 

K(s) = \fA^dr) iT+m > 


(19) 


( 20 ) 
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Expanding (flUl) in powers of W qq and averaging over DA, one obtains the 
effective one-loop partition function mm 

1 poo (jo 

Zuoop = — 2 l -(D 4 z) xx d 4 x e~ K ^ X (W qq , W qQ ), (21) 

where x is a a connected average of the product of two loops 


x = {W,s(A)Wqq(A)) - (W„(A))(W qq (A)). (22) 

The properties of x f° r different contour orientations of C qq and Cqq 
have been studied in m m ns!, and in [Q] it was found , that for the 
simplest case of the flat overlapping contours of opposite orientation one can 
approximate x as follows 

A ~ -j ^2 gx I>(-^a) = exp(—cr ren S) (23) 

where S'a is the area with subtracted area of loops qq, and a ren = cr is 
the string tension renormalized with account of sea quarks holes. 

We define the density of the sea quark holes in the confining film in 
Wqq(R, T), p = where AS is the area of the holes, in the case of zero 
m.f., and follow the development of p with the magnetic held. It is clear, that 
the increasing energy of the holes yields the increase of the effective string 
tension, which can be estimated from (ITT]) as 

Aa(eB) _ ^AE(B) _ AS^^B ^ 

a ~ P E(R ) S oR ’ ^ ' 

while the growth of p due to magnetic focusing in the case of B\\ should 
decrease effective string tension with Act pHeB). 

As a result one can write, taking into account, that 


Aa(eB = 0) 

CT 


AS Aa(eB) 

ri - Pi 

S CT 


Act (eB) — Aa(eB = 0) 

CT 


(25) 


where 


Act (eB) 

CT 


A S' 

— ( f± (eB)-MeB)), 


( 26 ) 
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(27) 


f±(eB) = 


e q B | 


<7 R 


and for f\\(eB) one has only the magnetic focusing effect, 


/||(efl) = £(efl)-l. (28) 

Note, that the signs of both terms (l27j) . (|28l) in Eq. (|26l) . are opposite. 
One must have in mind, that the term f\\(eB) is present for the parallel 
direction of the m.f., B = B||, while the second term on the r.h.s. in (1261) 
f±(eB) is active when magnetic field is perpendicular to the area. 


5 Perturbative gluon exchange in magnetic 
field 


We now turn to the gluon exchange interaction Voge in magnetic field, which 
was studied on the lattice in [5] [TB] [XT] and analytically in [IB], and exploited 
in [llJj to predict the meson mass behavior in m.f. 

ft was argued in [IS], that m.f. creates a screening effect in Voge due 
to the appearance of the quark loop contribution, which grows in m.f. in 
the same way, as the quark pair energy (ITU]) . This effect was known for a 
long time [20] and was exploited in [2Tj to predict the saturating effect in 
QED. Following this line in the framework of QCD in [TB] was obtained the 
one-loop Vqge with the dependence on m.f. in the form 


Vqge(Q) = 


167T 


a x 


Q 2 (1 + »sn f \e q B\ 

r f n Q 2 exp 




2\ eq B\ 


where T(z) = Q = q± + q \, and 


a 


(o) 


CM. a 


l + ^p\n Q2+Ml 


Mi 


/3ln 


— _• R = —N 

Q 2 + Mg ’ ^ 3 C ‘ 

A QCD 


(29) 


(30) 


When one is measuring Voge(R) on the lattice with R ~ (0.5 1) fin, 

one has gjj < a and small or vanishing qj_. Correspondingly one can expand 
T(z) and rewrite (j29l) as 


( V 0 GE(R))eB — (Voge(R)) 0. iy - , C — 9 ' (31) 

1 + Ca s bn a 
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and 


a 


(32) 


A a{eB) 
a s 


— Oi s + 


S 


1 + Ca s 


Ca s 
1 + Ca s 


Note the difference between the screening situation in QCD and QED. In 
QED there is no string, and hence no string direction R, and the exchange 
and the e + e _ , loops, transverse with respect to B, become heavy ( VeB, 
Eq. (fUJ) ) and this effect screens the Coulomb interaction in the transverse 
direction. 

In QCD the confining him (the string) defines the direction R, with the 
sea quark loop lying inside the him and hence one should have the screening 
effect as in (ITZll for R_LB and no screening in the case R||B, when sea quarks 
move in the loops along m.f. 

In this case, however, the focusing effect, £(eR) > 1, is acting, increasing 
the sea quark loop density p = =f- as — 1). This density is entering 

the general one-loop expression (130|) for cx s , where (3 = fio = ylYc ~ \ n f and 
the last two factors estimate the relative density of gluon and quark loops 
respectively. In our case the increased density of quark loops leads to the 
replacement in (1301) 


P > A) + A J3 0 — —A c — -rif — -n/(£(eB) — 1), (33) 

since only 1/3 the rif quark loops lies in the parallel to B position. 
Expanding in (l30li in powers of (£ — 1), one obtains 


Aaq 

cx s 


2 rif 

Wo 


(£-!)■ 


(34) 


One can see different signs of the m.f. action on a s in (l32l) and (J34]). 


6 Comparison to the effective action expan¬ 
sion 

We now turn to the general arguments, based on the expansion of the effective 
action S e ff, corresponding to (fTU|h namely we define as in [T7] . appendix D, 

^det(m q + D(A))} = exp (-S e f f ). 


8 


(35) 






The fourth order term in the expansion of S e ff in powers of constant field 
terms was obtained in [22] and generalized to the case of the superposition 
of magnetic held B and colorelecric held £ and colormagnetic 13 in [T7] . 

The 0(B 2 ) contribution has the form (see Eq. (D.5) from [T7]) 

sfif = + irSi + tr£l ~ ¥ e S' (36) 

Now taking into account, that the partition function Z (fl9l) is propor¬ 
tional to exp (—S e ff), one can immediately see, that in the case 5TB (i.e. 
B orthogonal to the Wilson loop surface, which was denoted above in the 
paper as the case of Bj_), one has 

exp (—S e ff) = exp(|const|(e.B) 2 5 2 ) > 1, (37) 

while in the case of Bn one obtains 

exp(— S e ff) = exp(—|const|(ei?) 2 5jj) < 1. (38) 

The string tension is obtained from the correlator of the colorelectric helds 

m 

o- = iy D(x)d 2 x , D(x,y) ~ (tr£i{x)^£i{y)</>). (39) 

Note, that exp(— S e ff) enters as a factor in the held averaging denoted 
by angular brackets (l39lh Hence Eqs. (1371) . (I38|) tell us, that 

&(£±, B) > cr(0, 0) (40) 

a(5| h B) <a(0,0) (41) 

in agreement with lattice measurements of [5]. 

On the lattice the relevant behavior for (5jj) and {£\) was found hrst in 
[22J for the SU( 2) group and in [17] for the real QCD and is in agreement 
with [5] and our results for Aa(B^) and Acr(T>j_) respectively. 

Note also, that in (l37jh (j38li the coefficient of 5 2 is 2.5 times bigger than 
that of £]_, which is qualitatively similar to our relations of Acr(B||) and 
A a(B±), following from /y and f±, Eqs. (12H1) and (1271) . 
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7 Comparison to the lattice data [5] 


To compare with numerical data one should fix the parameters, entering in 
our equations (l26]) - (j28j) . (132]) . (134|) . Actually, the relative density of the qq 
holes in the confinement area ^ is the only free parameter of our approach 
and we choose it as ^ = 0.15, i.e. we suggest that the holes of sea quark 
loops occupy ~ 15% of the whole area in absence of m.f. 

For f± (1271) and £(e£>), Eq. (TT6|) one should define the average value of 
(jj) in the lattice measurements, and we take it (j^) = 0 fm = 0.267 GeV, 
since a large part of measurement was done in the interval 0.5 fm < R < 1 fm. 
Correspondingly, er = 0.16 GeV 2 in (fTUlh and one obtains 


MsM = 0 . 31 a / —^— 2 . 
a VI GeV 2 ’ 


(42) 


where we have taken into account, that e q B = |eB for rif = 2 + 1. Now we 
turn to the function C,(eB) in (fT5lh which can be approximated as £(e5) = 
\J 1 + ^P-, and again with e q B = \eB, one has 


Acr(eS||) 

a 


As 

s 





( 43 ) 


The resulting curves of a (eB) for B± and B\\ are shown in Fig. Q] together 
with the lattice calculations of [5j (see Fig. 4 there at L — 40). One can 
see a quantitative agreement in both cases with our estimate ^ = 0.15, and 
important agreement can be seen in the low eB behavior, where (1431) yields 

quadratic growth <5cr(ell||) ~ ^ (f^) • 

We now turn to the case of a s (eB), Eqs. (T3T]l (t34|) . In (l3Tf one has for 
7if = 3, e q ~ e/2, C = 0.44 1 ^' y2 , and therefore with a s = 0.4, one has 


Aa s (eB 1 ) 


0.176(e5/GeV 2 


a, l + 0.176(eS/GeV 2 )’ 

and for a s (eB y), writing (134)) to all orders of 


Aa s (eB\\) 


Ota 


\ 

A+( 

' eB \ 
k 2k? ) 

2 

- 1 

14 '5 - \ 

A+ 

to 


(44) 


(45) 
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Figure 1: a(eB)/a(B = 0) for B± (upper red curve) and B\\ (lower red 
curve), in comparison with lattice data from [5] (green points) 

ft is clear that (145|) is valid for eB < 1 GeV 2 . 

The resulting curves of a(eB) for B± and B\\ are shown in Fig. [altogether 
with the lattice calculations of [1] (see Fig. 5 there at L = 40). One can 
find the same type of behavior, and again due to C,(eB), in our prediction for 
Aa s (eB\\) at small ( 73 ), which agrees well with the data of [5J. Our values 
for Aa s (eB ±) though are ~ 40% smaller in magnitude than the lattice data, 
but the general trend and sign are the same. Concluding, one can notice 
a qualitative agreement and, for the case of B\\, also a good quantitative 
agreement between our results and lattice measurements. 

8 Summary and conclusions 

Our discussion above is actually an attempt to qualitatively understand the 
dynamical mechanism beyond the a and a s dependence on m.f. We have 
identified two possible effects in the action of external magnetic held on 
confinement, which can act only through sea quarks loops. The first is the 
increasing production of loops in m.f. - the focusing effect. The second effect 
is the energy increase due to loop production, since they become effectively 
heavier in m.f., and this acts only when m.f. is perpendicular to the area 
surface. As a result one obtains different signs of combining effects; as shown 
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0 0.2 0.4 0.6 0.8 1 1.2 

eB, GeV 2 


Figure 2: a s (eB)/a s (B = 0) for B± (lower red curve) and B\\ (upper red 
curve), in comparison with lattice data from [5j (green points) 

in Fig. |T]and Fig. [2] and this corresponds to the lattice data [5j. 

To make quantitative comparison with lattice data of Eng m the only 
fitting parameter is p = which was taken as 0.15 . Note the difficulty 
in deriving it from the general theory eh Els], since the corresponding inte¬ 
grals are diverging and need regularization. The results for A a are in a fine 
agreement with [5j. 

We have calculated the screening of the a s (eB ) due to the quark pair 
creation in m.f., which occurs in Bi and has the same physical mechanism 
as in the qq energy growth due to m.f., Eq. (fill . The stimulated creation 
of the qq pairs in the case of B\\ leads to the increase of a s . The resulting 
forms of Aa s in Fig. [2] agree with lattice data [5j, and we have found the 
increase of a s in By due to the enhanced quark loop production, which is an 
antiscreening effect. 

The authors are grateful to Massimo D’Elia for discussions, suggestions 
and numerical data. The authors are grateful to M.A.Andreichikov, B.O.Kerbikov 
and members of the ITEP theory seminar for useful discussions. The finan¬ 
cial support of the RFBR grant 1402-00395 is gratefully acknowledged. 
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